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We investigate Bose-Einstein condensation for interacting bosons at zero and nonzero temper- 
ature. Functional renormalization provides us with a consistent method to compute the effect of 
fluctuations beyond the Bogoliubov approximation. For three-dimensional dilute gases, we find an 
upper bound on the scattering length a which is of the order of the microphysical scale - typically 
the range of the Van der Waals interaction. In contrast to fermions near the unitary bound, no 
strong interactions occur for bosons with approximately pointlike interactions, thus explaining the 
high quantitative reliability of perturbation theory for most quantities. For zero temperature we 
compute the quantum phase diagram for bosonic quasiparticles with a general dispersion relation, 
corresponding to an inverse microphysical propagator with terms linear and quadratic in the fre- 
quency. We compute the temperature dependence of the condensate and particle density n, and 
find for the critical temperature Tc a deviation from the free theory, ATc/T^ = 2.1a'n}^^ . For the 
sound velocity at zero temperature we find very good agreement with the Bogoliubov result, such 
that it may be used to determine the particle density accurately. 



I. INTRODUCTION 

For ultracold dilute non-relativistic bosons in three di- 
mensions, Bogoliubov theory gives a successful descrip- 
tion of most quantities of interest [l| . This approximation 
breaks down, however, near the critical temperature for 
the phase transition, as well as for the low temperature 
phase in lower dimensional systems, due to the impor- 
tance of fluctuations. One would therefore like to have 
a systematic extension beyond the Bogoliubov theory, 
which includes the fluctuation effects beyond the low- 
est order in a perturbative expansion in the scattering 
length. Such extensions have encountered obstacles in 
the form of infrared divergences in various expansions ■ 
Only recently, a satisfactory framework has been found 
to cure these problems [1, Ij] . 

Functional renormalization Q for the average action 
0) 9 systematically copes with the infrared problems 
by exploring the difficult long range behavior gradually 
by means of non-perturbative flow equations, which are 
based on an exact renormalization group equation 
For zero temperature they lead to a consistent description 
of Bose-Einstein condensation and the quantum phase 
transition for nonrelativistic bosons in arbitrary dimen- 
sion ^, .9]. A key ingredient is the generation of a term in 
the inverse propagator that is quadratic in the frequency 
~ vlJ^ . It characterizes the long time behavior of the 
full propagator at long distances, even if no quadratic 
frequency dependence is present in the microphysical or 
classical propagator. The coupling v is due to quan- 
tum fluctuations and appears during the flow from mi- 
crophysics to macrophysics. The extreme infrared limit 
is governed by a " relativistic" model where v dominates, 
with an enhanced (approximate) space-time symmetry 
SO{\^d) corresponding to Lorentz symmetry. In one or 
two dimensions the relativistic term ~ {i in the inverse 
propagator is crucial for a correct description of the low 
temperature behavior. In contrast, in three dimensions 



the flow towards the extreme infrared behavior is only 
logarithmic, such that the quantitative influence remains 
quite moderate for practical purposes. Nevertheless, the 
effective coupling v is needed in order to avoid the in- 
frared problems. 

In this paper, we extend this formalism to a nonva- 
nishing temperature. We present a quantitative rather 
accurate picture of Bose-Einstein condensation in three 
dimensions and find that the Bogoliubov approximation 
is indeed valid for most quantities. The same method 
can be applied for one or two dimensions, such that the 
present work can also serve as a test of the method. This 
is nontrivial, because insufficient truncations of the flow 
equations can lead to fake dependencies on the micro- 
scopic physics. 

For dilute non-relativistic bosons in three dimensions 
we find an upper bound on the scattering length a. This 
is similar to the "triviality bound" for the Higgs scalar 
in the standard model of elementary particle physics. As 
a consequence, the scattering length is at most of the or- 
der of the inverse effective ultraviolet cutoff A^^, which 
indicates the breakdown of the pointlike approximation 
for the interaction at short distances. Typically, A~^ is 
of the order of the range of the Van der Waals inter- 
action. For dilute gases, where the interparticle distance 
j^-i/3 jg larger than A~^ , we therefore always find a 

small concentration c = av}/^ . This provides for a small 
dimensionless parameter, and perturbation theory in c 
becomes rather accurate for most quantities. For typical 
experiments with ultracold bosonic alkali atoms one has 
A^i « 10"'^ cm, n^/^ sa 10'*cm"\ such that c < 10"^ is 
really quite small. 

Bosons with pointlike interactions can also be em- 
ployed for an effective description of many quantum 
phase transitions at zero temperature, or phase transi- 
tions at low temperature T. In this case, they correspond 
to quasi-particles, and their dispersion relation may dif- 
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fer from the one of non-relativistic bosons, uj = We 
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describe the quantum phase transitions for a general mi- 
croscopic dispersion relation, where the inverse classical 
propagator in momentum and frequency space takes the 
form Gq^ = —Suj — Vlu"^ (in units where the particle 
mass AI is set to 1/2). We present the quantum phase 
diagram at T = in dependence on the scattering length 
a and a dimensionless parameter v V/S'^, which mea- 
sures the relative strength of the term quadratic in lu in 
Gq^. In the limit S ^ (u — > oo) our model describes 
rclativistic bosons. 

A nonzero temperature T introduces a new scale into 
the problem. This modifies the effective scaling behav- 
ior. In particular, near the critical temperature T^, where 
the Bose-Einstein condensate dissolves, the long distance 
physics can be described by classical statistics. The non- 
relativistic bosons belong to the universality class of the 
three-dimensional 0(2) model, with the associated uni- 
versal critical exponents. The flow equations in the clas- 
sical regime are well studied and lead, for an appropri- 
ate level of the truncation, to very accurate results for 
the critical behavior [1, . Here we use a rather sim- 
ple truncation, but we obtain nevertheless a reasonable 
description of the non-analytical critical behavior. How- 
ever, the temperature range for the applicability of the 
universal critical behavior is found to be very small for 
the small values of c encountered for the dilute gases. 

The value of the critical temperature for interacting 
bosons cannot be computed within the Bogoliubov the- 
ory. We find a temperature shift ATc/Tc — 2.1an^/^, 
compared to the free theory. In contrast to other es- 
timates, which only take the classical fluctuations into 
account (only the zero Matsubara frequency), we include 
the full quantum statistics (all Matsubara frequencies). 
We finally compute the sound velocity at T = 0. This 
quantity agrees to high precision with the Bogoliubov re- 
sult. For a known scattering length a, it can therefore be 
used as a precise measure of the density. 

Our paper is organized as follows. As a starting point, 
we specify our microscopic model in Sec. |TI1 In the fol- 
lowing Sec. mil we recall the method of functional renor- 
malization, explain our truncation of the effective aver- 
age action and also describe the projection of the exact 
flow equation onto that truncation. The flow equations 
in the vacuum are investigated in Sec. IIVI and an up- 
per bound for the scattering length is derived. In Sec. 
IVl we explain our method to determine the density and 
also study the quantum phase diagram at zero tempera- 
ture. The effects of a non-vanishing temperature on the 
condensate are taken into account in Sec. IVIi while Sec. 
I VIII gives our result for the critical temperature in depen- 
dence of the interaction strength. Finally, we investigate 
the sound velocity and draw conclusions in Sees. IVIIII 
andim 

The Appendices |X] and |B] contain a motivation of our 
truncation in terms of a systematic derivative expansion 
and an analysis of symmetries constraining the form of 
the effective action. Our explicit results for the flow equa- 
tions of the effective potential and the kinetic coefficients 



are shown in Appendices [Cl and [P] Appendix [E] analyzes 
the dispersion relation. 

II. MICROSCOPIC MODEL 

Our microscopic action describes nonrelativistic 
bosons, with an effective interaction between two par- 
ticles given by a contact potential. It is assumed to be 
valid on length scales where the microscopic details of the 
interaction are irrelevant and the scattering length is suf- 
ficient to characterize the interaction. The microscopic 
action reads 

^[0]=^ (^5.-1^9,2- A-a)0-HiA(0* 0)2}, (1) 
with 

x= {t,x), J ^ Jj dr J Sx. (2) 

The integration goes over the whole space as well as over 
the imaginary time r, which at finite temperature is in- 
tegrated on a circle of circumference /? = l/T accord- 
ing to the Matsubara formalism. We use natural units 
h = ks = 1. We also scale time and energy units with 
appropriate powers of 2M, with M the particle mass. In 
other words, our time units are set such that effectively 
2M =1. In particular, we use the rescaled chemical 
potential a — 2Mfi and T stands for the temperature 
multiplied by 2M. In these units time has the dimension 
of length squared. For standard non-relativistic bosons 
one has = and = 1 , but we also consider quasipar- 
ticles with a more general dispersion relation described 
by nonzero V. 

After Fourier transformation, the kinetic term reads 

U*{q){lSq^ + Vql+f)m. (3) 

Jq 

with 

At nonzero temperature, the frequency go ~ =^ 2-KTn 
is discrete, with 




while at zero temperature this becomes 

The dispersion relation encoded in eq. ^ obtains by 
analytic continuation 

Su + Vlo'^ = f/2M. (7) 
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In this paper, we consider homogeneous situations, i.e. 
an infinitely large volume without a trapping potential. 
Many of our results can be translated to the inhomo- 
geneous case in the framework of the local density ap- 
proximation. One assumes that the length scale rele- 
vant for the quantum and statistical fluctuations is much 
smaller than the characteristic length scale of the trap. 
In this case, our results can be transferred by taking 
the chemical potential position dependent in the form 
a (x) = 2M{fj, — Vt{x)), where Vt{x) is the trapping po- 
tential. 

The microscopic action ([1} is invariant under the global 
U{1) symmetry which is associated to the conserved par- 
ticle number, 

</.^e"(/.. (8) 

On the classical level, this symmetry is broken sponta- 
neously when the chemical potential a is positive. In this 
case, the minimum of —a(j)*(j) + ^X{(j>*(j))^ is situated at 
= J. The ground state of the system is then charac- 
terized by a macroscopic field 00 j with 0q0o = /Oq = j- 
It singles out a direction in the complex plane and thus 
breaks the U{1) symmetry. Nevertheless, the action it- 
self and all modifications due to quantum and statistical 
fluctuations respect the symmetry. For V = and = 1, 
the situation is similar for Galilean invariance. At zero 
temperature, we can perform an analytic continuation to 
real time and the microscopic action ([T]) is then invari- 
ant under transformations that correspond to a change 
of the reference frame in the sense of a Galilean boost. 
It is easy to see that in the phase with spontaneous U{1) 
symmetry breaking also the Galilean symmetry is bro- 
ken spontaneously: A condensate wave function, that is 
homogeneous in space and time, would be represented in 
momentum space by 



(9) 



Under a Galilean boost transformation with a boost ve- 
locity 2q, this would transform according to 

^ (j){uj ~ f,p- q) 

^q^oi2n)U(^Hp~-q)Siu;^f). (10) 

This shows that the ground state is not invariant under 
such a change of reference frame. This situation is in 
contrast to the case of a relativistic Bose-Einstein con- 
densate, like the Higgs boson field after electroweak sym- 
metry breaking. A relativistic scalar transforms under 
Lorentz boost transformations according to 

<j>ipn^Hi^-'npn, (n) 

such that a condensate wave function 

^0 (2^)^ <5(4) {pn 00 {2nr <S('^ ((A"!)^, p^) 

-0o(27r)4^W(p^) (12) 

transforms into itself. We will investigate the implica- 
tions of Galilean symmetry for the form of the effective 



action in app. [X] An analysis of general coordinate in- 
variance in nonrelativistic field theory can be found in 

El. 



III. NON-PERTURBATIVE FLOW EQUATIONS 

A. Functional Renormalization Group and Flow 
equation 

We start with a functional integral representation of 
the grand canonical partition function 



Z = Tr e 



-l3iH-f,N) 



Dxe 



~S[x] 



(13) 



In this paper, we work with the formalism of quantum 
statistics for many particle problems. In contrast to clas- 
sical statistics, the fields %(''": x) are parametrized by an 
Euclidean time variable t in addition to the space vari- 
able X. This Euclidean time is wrapped up on a circle of 
circumference /3 = y, such that the fields x{t,x) live on 
a (generalized) torus. The microscopic action S[x\ is of 
a form similar to ([T]). 

We generalize eq. (|13p by introducing a source J for 
the fields x write 



(14) 



The rt-point correlation functions can now be obtained 
by functional differentiation of Z[J], while W[J] gener- 
ates the connected n-point functions. For example, the 
expectation value of x obtains from 



(x) = {x{x)) = 



6W[J] 
5J{x) ' 



(15) 



The thermodynamic potential '^q associated with the 
grand canonical partition function is given by 



1. 



-\nZ 



1 



W[J = 0]. 



The effective action r[0] is defined as 



j(j) 



J=.h: 



where Jcx is obtained by the inversion of 
6W[J] , 



SJ 



(16) 



(17) 



(18) 



It is straightforward to show 



sr\, 



J. In the absence of 



a source J, we obtain for the thermodynamic potential 
of the grand canonical partition function 



r[0cq] 



(19) 
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where (/)oq is defined by 



= 0. 



(20) 



The effective action r[0] is the generating functional of 
the one-point irreducible correlation functions and in a 
sense, its precise knowledge corresponds to the solution 
of the theory. 

Our method determines r[0] with the help of an exact 
flow equation. For that purpose, we include an infrared 
cutoff term ASk [x] in eq. and define 



^W,[J] ^ J £)^g-S[x]-ASfc[x]+/ Jx_ 

In Fourier space, the cutoff term reads 



and has the properties 



Rkiq)x*{q)x{q) 



(21) 



(22) 



The effective average action is defined as a modified Leg- 
endre transform of Wk 



-Wk[J] 



J(t) 



ASk[<t>\ 



(24) 



It has the important property, that it interpolates be- 
tween the microscopic action S[(j)\ and the full effective 
action r[(/)] 



TkWi -^S[4>] (fc^oo), 

Tk[4>] ->r[0] (fc->o). 



(25) 



In physical terms Tk [4>\ is the effective action in the pres- 
ence of an infrared cutoff at a momentum scale k. Only 
fluctuations with momenta larger than k are included. 
For example, a finite volume ~ -p- of the system under 
consideration would lead to an situation that is described 
by Ffe [(/)]. 

Our method to determine Vk[<l>\ (and for k ^ also 
relies on the existence of an exact flow equation 



SfeFfc 



Tr(ri^) 



Rk) ^dkRk 



(26) 



Here the trace operation includes a momentum integra- 
tion J^, as well as a sum over internal indices i = 1,2, 
according to the two real components in the decomposi- 
tion (j){x) = -i=(0i(a;) +i(f)2{x)). On the r.h.s. of (gSl), 

(2) w. 



T)^' is the second functional derivative of Tk[4>] 



^{-q) 



(27) 





oo ( k - 


-y oo), 


truncate 


Rk{q) 


=^ fc2 {q- 






Rk{q) 


-^0 (fc- 


-0). 


(23) 



It is therefore a matrix in internal and momentum space. 
Correspondingly, Rk in eq. (pS)) stands for Rk{q)6ij6{q — 
p). The flow equation ([26]) describes the evolution of the 
effective average action with the cutoff scale k. 

The functional differential equation ((26)l is hard to be 
solved exactly. In principle, Ffc is described by infinitely 
many couplings. We will use here an approximation with 
only a finite number of couplings. This is achieved by an 
ansatz for a specific form of Tk[(t)]- In the absence of 
anomalies the effective action r[0] is invariant under the 
same symmetries as the microscopic action S[4>\. This 
holds also for the effective average action F^ , provided 
that the cutoff term AS'fc[(/)] is also invariant. Our ansatz 
will respect all symmetries of the classical action. 



B. Truncation 

Approximate solutions of the exact ffow equations ob- 
tain from a truncation of the general form of the effective 
action. We use here terms with up to two derivatives and 



+2V{a -ao)4>* {Or - A) + U{p, a) } , (28) 



with p — (t>*(t>- This particular form is motivated by a 
more systematic derivative expansion and an analysis of 
symmetry constraints (Ward identities) in appendix [Xl 
We introduce the renormalized fields </> = A-^^^cj), p — 
Ap, the renormalized kinetic coefficients S = V = ^ 
and we express the effective potential in terms of the 
renormalized invariant p, with 

Uip,a) = U{p,a). 

This yields 

+2Via ~ao)4>* [dr - A) + U{p, a) ) . (30) 



(29) 



For the effective potential, we use an expansion around 
the fc-dependent minimum po{k) of the effective potential 
and the fc-independent value of the chemical potential ctq 
that corresponds to the physical particle number density 
n. We determine po{k) and ctq by the requirements 

{dpU){po{k),ao) =0 for aUfc 
~{d„U)(po,<JQ) atfc^O. (31) 

More explicitly we take a truncation for U{p, a) of the 
form 

U{p,a) = U{po,ao) - nk{(T - ao) 



CTo)) (p - Po) 



+ -{X + (3{a-ao))ip-pof. (32) 
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In the symmetric phase we have po = 0, while in the 
phase with spontaneous symmetry breaking, we have 
rri^ — 0. In summary, the flow of for fixed a — aQ 
is described by four running renormalized couplings po, 
A, S and V. In addition, we need the anomalous di- 
mension 7] = —kdklnA. A computation of n requires a 
flow equation of rife, which involves the couplings linear 
in cr — (To 5 namely a and /3. The pressure is calculated 
by following the fc-dependence of the height of the mini- 
mum pk = -U{po,ao). All couplings po, \ S, V, A, rifc, 
Pk, a, (3 depend on k and T. The physical renormalized 
couplings obtain for fc — > 0. They specify the thermody- 
namic potential C/(pojfo) as well as suitable derivatives 
of the potential and the correlation function. 

The "initial values" at the scale A; = A are determined 
by the requirement 



(33) 



functional derivative of Tk- For a real constant back- 
ground field 4'i{x) = \/2p, 4'2{x) = the latter becomes 
diagonal in momentum space 



For our truncation one has at a = (Tq 



(39) 



G- 



A 



f + Vq^ + U' + 2pU" 



Sqo 



-Sqo 
Vq'o^ 



U' 



(40) 

Here, primes denote derivatives with respect to p (not 
p). In the phase with spontaneous symmetry breaking, 
the infrared cutoff in the flow equation ([26|) adds to the 
diagonal term in (|^ a piece A{k'^ — (f)6{k'^ — if). This 
effectively replaces (7^ — > fc^ in eq. (|26p whenever < k^, 
thus providing for an efhcient infrared regularization. 



using the microscopic action S[(j)\ in eq. ([1]). This implies 
the initial values 



'0'o)cto, 



aA = -l, /3a = 0. (34) 



We remain with the free microscopic couplings Xa, Sa = 
Sa, Va — Va- The coupling Aa will be replaced by the 
scattering length a in the next section. We further choose 
units for r where S'a = 1. Then our second free coupling 



02 



^aA^ 



In consequence, besides the thermodynamic variables T 
and fJo our model is characterized by two free parameters, 
a and v. Often, we will concentrate on "standard" non- 
relativistic bosons with a linear r derivative, such that 
V = 0. The scattering length a remains then the only 
free parameter. In the vacuum, where T ^ n ~ 0, this 
sets the relevant unit of length. 

Finally, we have for the infrared cutoff 



C. Non-perturbative flow equations 

We project the flow equation of the effective average 
action onto equations for the coupling constants by us- 
ing appropriate background fields and taking functional 
derivatives. The flow equation for the effective potential 
obtains by using a space- and time-independent back- 
ground field in eq. ([^ . with t = ln(A:/A) 



dtU\^ = kdkU\^ = C = - ivGkdt{Ark). (41) 

"'9 



(35) The propagator Gk is here determined from 



Gfc ' = G-^ 



with 



Pii 
P21 
P22 



A 



Vql 



Tk 
Tk 



^A 



P21, 



^12 
P22 



(42) 



-Pt 



e + Vql 



U' + 2pU" + 2V{a-ao)f, 
Sqo + 2V{a - (Jo)qo, 

U' + 2V{a- ao)f. (43) 



A5fc= / Arrk{-^)4>^ / </>*rfc(-A),/) 



We choose the optimized cutoff function [ij] 

rk{p^) = (fc2 - p2 _ ^2)^^^2 _ p2 _ ^2)^ 



(36) 



(37) 



where we recall, that m? = in the regime with spon- 
taneous symmetry breaking. It is convenient to work 
with real fields 0i,2(a^), (t'ix) = '^{4'i{^) +'i4'2{x)), with 
Fourier components 



Again primes denote a differentiation with respect to p. 
We switch to renormalized fields by making a change of 
variables in the differential equation (|4T|) 



ripU'. 



(44) 



We can now derive the flow equations for the couplings 
pti{k) and X{k) by appropriate differentiation of ((44|) with 
respect to p. The flow equation for U is given more ex- 
plicitly in appendix [Cl Differentiation with respect to a 
yields the flow of n^, a, (3. We use in detail 



0,(T,f)= / e'«-0,(9) 



e'^^°^+^"^j{qo,q). (38) 



The inverse propagator for the fields (p becomes a 2 x 2 
matrix in the space of 4>i and 4>2, given by the second 



dt 



= ^^{dlU){po,ao) 

= {dptU)ipQ,<Jo) + idp)ipo,<Jo)dtPo, 
= 9lC\ 



Pa,iyo 



2r,A, 



(45) 
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and similar for the flow of a and (3, 
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FIG. 1: (Color online) Flow of the interaction strength A with 
the scale parameter k for different start values, corresponding 
to Aa = 10^ (dotted), Aa = 0.044 (solid) and Aa = 0.0026 
(dashed). The first case is plotted for zero density (n — 0) 
only, while the last two cases are plotted also for unit density 
{n — 1). The curves n = and n — 1 are identical within the 
plot resolution. The solid and the dashed curve correspond 
to a = 10""^ and a = 10~^, respectively. 
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FIG. 2: (Color online) Flow of the minimum of the effective 
potential for n = 1. The parameters for the solid and the 
dashed curves are the same as in figure [T] 



where we recall, that d^U = in our truncation. To 
determine the flow equation of po, we use the condition 
that U'{po) = for all k, and therefore 



iidpU)ipo,ao)^0, 
{dpdtU){pa, fJo) + (52C/)(po, ctq) dtpa 0, 



.(46) 



We show the flow of A and po in figs. [UOfor n = 1, T = 
and different values of Aa (with -5 = 0). The change in 
Pq is rather modest. This will be different for nonzero 
temperature. 

The flow of rik is given by 



dt 



d 



{-d„U){pQ,aQ) 



-{d„dtU){po,ao) - {dpdaU){po,(Jo) dtpo 



= -daC 



adtPo, 



(47) 



dt 



dt 



d 
di 

{dpdadtU){po,ao) + (9p9a?7)(po, cto) 3tPo 

d 

di 



{dpdaU){po,ao) 



{dld,U){po,<Jo) 



{dldadtU){po,(To) + {dpdaU){po,ao)dtpo, 



277/3, 



(48) 



where the last equation holds since dpdaU = in our 
truncation. 

For a derivation of 77 = —{dtA)/A and the flow equa- 
tions for S and V, we have to evaluate the flow equation 
(PS)) for a background field depending on go and q. We 
use an analytic continuation go = 'i-'^ and obtain the flow 
equation for S from 



dtiSA) = -in 



-1 d 



duj (-w,0) (5(/'i(a;,0) 



dtTk 



j=0 



with four-volume O = ^ J^. The projection prescription 
for V is 



dtivA)^-n 

and similar for A 
d 



1 d 



dtTk 



a)=0 



(50) 



(51) 



After the functional differentiation, we evaluate the ex- 
pressions (|49| . ([50|) . ([5T|) at homogeneous background 
fields. These calculations are a little intricate, but stan- 
dard and straightforward in principle. More explicit flow 
equations are given in app. [D] Eventually, it is always 
possible to perform the Matsubara sums and also the 
spatial momentum integration analytically. In figure [3] 
we show the flow of A, S and V at zero temperature 
and for density n = 1. The kinetic coefficient A starts 
on the large scale with A = 1, increases a little around 
k = v}!'^ and saturates to a constant. In contrast, the 
coefficient S starts to decrease after a short period of in- 
crease with A. For very tiny scales fc, S would finally go 
to zero. The frequency dependence of the propagator is 
then governed by the quadratic frequency coefficient V . 
In three spatial dimensions, however, this decrease of S 
is so slow that it is not relevant on the length scales of 
experiments. This is one of the reasons why Bogoliubov 
theory, which neglects the appearance of V , describes ex- 
periments with ultracold bosonic quantum gases in three 
dimensions with so much success. The coefficient V is al- 
ways generated in the phase with spontaneous symmetry 
breaking 0]. Its fc-dependence is also shown in figure [31 
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FIG. 3: (Color online) Flow of the kinetic coefficients A 
(solid), S (dashed) and V for a scattering length a = 10"'^, 
temperature T = and density n = 1. 



IV. UPPER BOUND FOR THE SCATTERING 
LENGTH 

The vacuum is defined to have zero temperature T = 
and vanishing density n = 0, which also implies — Q. 
The interaction strength A at the scale fc = determines 
the four point vertex at zero momentum. It is directly 
related to the scattering length a for the scattering of two 
particles in vacuum, which is experimentally observable. 
We therefore want to replace the microscopic coupling Aa 
by the renormalized coupling a. In our units (2Af = 1), 
one has the relation 



a = — A(fc = 0,r = 0,n = 0). 
Stt 



(52) 



The vacuum properties can be computed by taking for 
T ~ Q the limit n — > 0. We may also perform an equiva- 
lent and technically simple computation in the symmetric 
phase by choosing m?{k = A) such, that w?{k ^ 0) = 0. 
This guarantees that the boson field </() is a gap-less prop- 
agating degree of freedom. 

We first investigate the model with a linear r- 
derivative, S'a = 1, Va = 0. Projecting the flow equation 
([26)1 . we find the following equations: 



dtrn^ = 



A2\ 



^2)3/2 



Qik^-m^). (53) 
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FIG. 4: (Color online) Scattering length a in dependence on 
the microscopic interaction strength Aa (solid). The asymp- 
totic maximum Omax = |x is also shown (dashed). 



The propagator is not renormalized, dtS — dtV — dtA — 
0, 77 = 0, dtOi = 0, and one finds dtUk = 0. The coupling 
/? is running according to 



dtf} = 



1 



1 



3 



fc4 7r2S' 



(54) 



Since f3 appears only in its own flow equation, it is of no 
further relevance in the vacuum. Also, no coupling V is 
generated by the flow and we have therefore set V — 
on the r.h.s. of eqs. ([55)1 and ([51]) . 

Inserting in eq. (j53p the vacuum values = and 

= 1, we find 



67r2 



rA^ 



The solution 



A(fc) = 



(55) 



(56) 



tends to a constant for fc ^ 0, Aq = A(fc = 0). The 
dimensionless variable A = ^ goes to zero, when k goes 
to zero. This shows the infrared freedom of the theory. 
For fixed ultraviolet cutoff, the scattering length 



1 



Ao ^ 

Aa 37T 



(57) 



as a function of the initial value Aa, has an asymptotic 



maximum 



Stt 
4A' 



(58) 



The relation between a and Aa is shown in fig. |4l 

As a consequence of eq. ([58| , the nonrelativistic bosons 
in d = 3 are a "trivial theory" in the sense that the 
bosons become noninteracting in the limit A — > 00, where 
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a ^ 0. The upper bound ((58)) has important practical 
consequences. It tells us, that whenever the "macrophys- 
ical length scales" are substantially larger than the mi- 
croscopic length A~^, we deal with a weakly interacting 
theory. As an example, consider a boson gas with a typ- 
ical inter-particle distance substantially larger than A^^. 
(For atom gases A~^may be associated with the range 
of the Van der Waals force.) We may set the units in 
terms of the particle density n, n = 1. In these units A 
is large, say A = 10^. This implies a very weak interac- 
tion, a < 2.5 • lO^''. In other words, the scattering length 
cannot be much larger than the microscopic length A~^. 
For such systems, perturbation theory will be valid in 
many circumstances. We will find that the Bogoliubov 
theory indeed gives a reliable account of many proper- 
ties. Even for an arbitrary large microphysical coupling 
(Aa — > oo), the renormalized physical scattering length a 
remains finite. 

Let us mention, however, that the weak interaction 
strength does not guarantee the validity of perturbation 
theory in all circumstances. For example, near the criti- 
cal temperature of the phase transition between the su- 
perfluid and the normal state, the running of X(k) will 
be different from the vacuum. As a consequence, the 
coupling will vanish proportional to the inverse correla- 
tion length as T approaches Tc, A ~ T~^^~^. In- 
deed, the phase transition will be characterized by the 
non-perturbative critical exponents of the Wilson-Fisher 
fixed point. Also for lower dimensional systems, the up- 
per bound (|58p for Aq is no longer valid - for example the 
running of A is logarithmic for d — 2. For our models 
with Va ^ 0, the upper bound becomes dependent on 
Va- It increases for Va > 0. In the limit Sa 0, it 
is replaced by the well known "triviality bound" of the 
four dimensional relativistic model, which depends only 
logarithmically on A. Finally, for superfluid liquids, as 
''He, one has n ^ A^, such that for a ~ A~^ one finds a 
large concentration c. 

The situation for dilute bosons seems to contrast with 
ultracold fermion gases in the unitary limit of a Feshbach 
resonance, where a diverges. One may also think about a 
Feshbach resonance for bosonic atoms, where one would 
expect a large scattering length for a tuning close to res- 
onance. In this case, however, the effective action does 
not remain local. It is best described by the exchange of 
molecules. The scale of nonlocality is then given by the 
gap for the molecules, uim- Only for momenta < m\j 
the effective action becomes approximately local, such 
that A = vriM for our approximation. Close to resonance, 
the effective cutoff is low and again in the vicinity of a^^. 



V. QUANTUM PHASE DIAGRAM 

A. Different methods to determine the density 

The density sets a crucial scale for our problem. Its 
precise determination is mandatory for quantitative pre- 



cision. We will discuss two different methods for its de- 
termination and show that the results agree within our 
precision. For T = 0, we also find agreement with the 
Ward identity n — pQ. 

The first method is to derive flow equations for the den- 
sity. This has the advantage that the occupation numbers 
for a given momentum p are mainly sensitive to running 
couplings with fc^ 
ism, the density is defined by 



p^. In the grand canonical formal- 



We can formally define a fc-dependent density rifc by 



dan 



(59) 



-(a,C/)(po,fTo). (60) 



The flow equation for is given in eq. (|47|) and the 
physical density obtains for A: = 0. The term 9o-C|p ^ 
that enters eq. ([47l) is the derivative of the flow equation 
(|¥T|) for U with respect to cr. To compute it, we need the 
cr-dependence of the propagator Gk in the vicinity of gq. 
Within a systematic derivative expansion, we use the ex- 
pansion of U{p,a) and the kinetic coefficients Zi and Z2 
to linear order in (cr — 170)7 ^ described in appendix [K\ 
Here, Zi(p,a) and Zii^p^u) are the coefficient functions 
of the terms linear in the r-derivative and linear in A, 
respectively. No reasonable qualitative behavior is found, 
if the linear dependence of Z\ and Zi on (cr — ctq) is ne- 
glected. Also, the scale dependence of a and /3 are quite 
important. The flow equations for a and /3 can be ob- 
tained directly by differentiating the flow equation of the 
effective potential with respect to a and p, cf. eq. (|48p . 
The situation is more complicated for the kinetic coef- 



ficients z'f '^ = daZi{po,aQ) and Z^"' = daZ2{po, ctq). 
Their flow equations have to be determined by taking 
the cr-derivative of the flow equation for Zi(p, cr) and 
Z2{p,a). As discussed in app. El we use in this paper 

the approximation z'f'^ = z'^^ = 2V = 2Vi{po, ctq). 

As a check of both our method and our numerics, we 
also use another way to determine the particle density. 
This second method is more robust with respect to short- 
comings of the truncation, but less adequate for high pre- 
cision calculations as needed e.g. to determine the con- 
densate depiction. The second method determines the 
pressure p = —U{po,ao) as a function of the chemical 
potential ctq. Here, the effective potential is normalized 
by U{po = 0, cro) = at T = 0, n = 0. The flow of the 
pressure can be read of directly from the flow equation 
of the effective potential and is independent of the cou- 
plings a and f3. We calculate the pressure as a function 
of cr and determine the density n = -^p by taking the cr- 
derivative numerically. It turns out that p is in very good 
approximation given by p = ccr^, where the constant c 
can be determined from a numerical fit. The density is 
thus linear in cr. 

At zero temperature and for v = Q, we can additionally 
use the Ward identities connected to Galilean symmetry. 
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FIG. 5: (Color online) Pressure and density as a function 
of the chemical potential at T = 0. We use three different 
methods: n = —dcrUmin from the flow equation (triangles), 
n = po as implied by Galilean symmetry (stars) and n — daP, 
where the pressure p — —U (boxes) was obtained from the 
flow equation and phenomenologically fitted by p = 56.5(t^ 
(solid lines). Units are arbitrary and we use a = 3.4 • 10~^, 
A = fO^ 



which yield n = pq. We compare our methods in figure 
[5] and find that they give numerically the same result. 
We stress again the importance of a reliable method to 
determine the density, since we often rescale variables by 
powers of the density to obtain dimensionless variables. 



B. Condensate and depletion density for T = 

We want to split the density into a condensate part 
ric and a density for uncondensed particles or " depletion 
density" nd ~ n — ric- For our model the condensate 
density is given by the "bare" order parameter 



= Po = Po(fc = 0). 



(61) 



In order to show this, we introduce occupation numbers 
n(p) for the modes with momentum p with normalization 



n{p} 



(62) 



One formally introduces a p dependent chemical potential 
a{p) in the grand canonical partition function 



g-i .„„l^l ^ Tre 

with three dimensional volume Q3 
define the occupation numbers by 



-fj{H-ihLa{p)n{p)) 



(63) 



n{p} 



6 



1 



Scrip) 



L. Then one can 

Jx 

(64) 



This construction allows us to use fc-dependent occupa- 
tion numbers by the definition 



nk{p) 



1 



5cj{p) pn^ 



tt)=4)a,(j(p)=ao 



(65) 



One can derive a flow equation for this occupation num- 
ber Ukip) [il: 



dknk{p) 



1 



l_t 

2 (5(7(p) /?r!3 

d 5 1 



Tr 



r[0,(7] 



00, fo 



(66) 



We split the density occupation number into a 5- 
distribution like part and a depletion part, which is reg- 
ular in the limit p ^ 



nkip) = nc.kS{p) + nd,k{p)- 



(67) 



One can see from the flow equation for nk(j)} that the 
only contribution to dkUc^k comes from the second term 
in equation ((66|) . Within a more detailed analysis [l^ 
one finds 



dkrick = dkpo.k- 



(68) 



We therefore identify the condensate density with the 
bare order parameter 



nc = Po ^ 



Po 
A 



(69) 



Correspondingly, we define the fc-dependent quantities 



nc,k = Po,k, nk = ric^k + nd,k 



(70) 



and compute rid = nd{k = 0) by a solution of its flow 
equation. 

Even at zero temperature, the repulsive interaction 
connected with a positive scattering length a causes a 
portion of the particle density to be outside the con- 
densate. From dimensional reasons, it is clear, that 
nd/n = (n — n^jn should be a function of an^/"^. The 
prediction of Bogoliubov theory or, equivalently, mean 
field theory, is ndjn — ^^(an^/'^)^/^. We may deter- 
mine the condensate depletion from the solution to the 
flow equation for the particle density, n — nk=o, and 
nc = Po = Po{k = 0). 

From Galilean invariance for T = Q and 5 = 0, it 
follows that 



n 



Po 



PO 



A' 



(71) 



with A = A{k = 0). This gives an independent deter- 
mination of ric- In flgure[6]we plot the depletion density 
obtained from the flow of n and pq over several orders of 
magnitude. Apart from some numerical fluctuations for 
small an^/^, we flnd that our result is in full agreement 
with the Bogoliubov prediction. 



C. Quantum phase transition 

For T — SL quantum phase transition separates the 
phases with po = and pq > 0. In this section, we in- 
vestigate the phase diagram at zero temperature in the 
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FIG. 6: (Color online) Condensate depletion (n — nc)/n as a 
function of the dimensionless scattering length an}^^ . For the 
solid curve, we vary a with fixed n — 1, for the dashed curve 
we vary the density at fixed a = 10~*. The dotted line is the 
Bogoliubov-Result {n — nc)/n = j^{an^^^)^^^ for reference. 
We find perfect agreement of the three determinations. The 
fluctuations in the solid curve for small an}'^^ are due to nu- 
merical uncertainties. Their size demonstrates our numerical 
precision. 



cube spanned by the dimensionless parameters tr = -p-, 
a — aA and v — -^A^. This goes beyond the usual 

phase transition for nonrelativistic bosons, since we also 
include a microscopic second r-derivative ~ v, and there- 
fore models vifith a generalized microscopic dispersion re- 
lation. For non-vanishing v (i.e. for a nonzero initial 
value of Vi with V2 = Vij = in app. E|, the GaHlean 
invariance at zero temperature is broken explicitly. For 
large v, we expect a crossover to the " relativistic" 0(2) 
model. If we send the initial value of the coefficient of 
the linear r-derivative Sa to zero, we obtain the limit- 
ing case V 00. The symmetries of the model arc now 
the same as those of the relativistic 0(2) model in four 
dimensions. The space-time- rotations or Lorentz sym- 
metry replace Galilean symmetry. 

It is interesting to study the crossover between the two 
cases. Since our cutoff explicitly breaks Lorentz symme- 
try, we investigate in this paper only the regime v ^ 1. 
Detailed investigations of the flow e quat ions for v 00 
can be found in the literature [1, [lol. Il4l| . The phase di- 
agram in the a — v plane with a = 1 is shown in figure 
[71 The critical chemical potential first increases linearly 
with V and then saturates to a constant. The slope in 
the linear regime as well as the saturation value depend 
linearly on a for a < 1. 

At r = 0, the critical exponents are everywhere the 
mean field ones {rj = 0, v = 1/2). This is expected: It 
is the case for w = and for V — 00 the theory 

is equivalent to a relativistic 0(2) model in d = 3 -f 1 
dimensions. This is just the upper critical dimension of 
the Wilson-Fisher fixed point . 

From section 3 we know that for v — the parameter 
a is limited to a < ^ ~ 2.356. For v — and a small 
scattering length a — > 0, a second order quantum phase 
transition divides the phases without spontaneous sym- 
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FIG. 7: Quantum phase diagram in the a-v plane for a — 1. 
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FIG. 8: (Color online) Quantum phase diagram in the a-a 
plane for v = 1 (dotted), v = 0.01 (dashed) and v = (solid). 



metry breaking for a < from the phase with a finite 
order parameter po > for cr > 0. It is an interest- 
ing question, whether this quantum phase transition at 
a — 0,v — also occurs for larger scattering length a. 
Wc find in our truncation that this is indeed the case for 
a large range of a, but not for d > 1.55. Here, the criti- 
cal chemical potential suddenly increases to large positive 
values as shown in fig. [S] For v > this increase happens 
even earlier. (For a truncation with Vi = 0, the phase 
transition would always occur at a — 0.) We plot the 
a — d plane of the phase diagram for different values of 
V in figure [5] The form of the critical line can be under- 
stood by considering the limits v ^ as well as d ^ 0. 

For a fixed chemical potential, the order parameter po 
as a function of a goes to zero at a critical value ac as 
shown in fig. [HI This happens in a continuous way and 
the phase transition is therefore of second order. For 
(T 0, we find Oc — 1.55A^^. We emphasize, however, 
that ttc is of the order of the microscopic distance A~^. 
Universality may not be realized for such values, and the 
true phase transition may depend on the microphysics. 
For example, beyond a critical value for the repulsive 
interaction, the system may form a solid. Ultracold atom 
gases correspond to metastable states which may lose 
their relevance for a — > A~^. For v > and cr ^ A^ the 
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FIG. 9: Quantum phase transition for fixed chemical potential 
(7 = 1, with A — 10^. The density po = n as a function of 
the scattering length a goes to zero at a critical OcA = 1.55, 
indicating a second order quantum phase transition at that 
point. 



phase transition occurs for OcA <C 1 such that universal 
behavior is expected. 
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FIG. 10: (Color online) Density nja^^'^ (solid) and order pa- 
rameter po/a^^^ (dashed) as a function of the temperature 
T/a. The units are arbitray with a = 2 • 10"* and A = 10^. 
The plot covers only the superfluid phase. For higher tem- 
peratures, the density is given by the thermal contribution 
n — riT only. 



VI. TEMPERATURE DEPENDENCE OF 
CONDENSATE 

So far, we have only discussed the vacuum and the 
dense system at zero temperature. A non vanishing tem- 
perature T will introduce an additional scale in our prob- 
lem. For small T <^ r?!'^ we expect only small correc- 
tions. However, as T increases the superfluid order will 
be destroyed. Near the phase transition for T ^ Tc and 
for the disordered phase for T > Tc, the characteristic 
behavior of the boson gas will be very different from the 
T ^ limit. 

For T > the particle density n receives a contribution 
from a thermal gas of bosonic (quasi-) particles. It is no 
longer uniquely determined by the superfluid density po- 
We may write 

n = po + UT (72) 

and observe, that nx = is enforced by Galilean sym- 
metry only for T = 0, Va = 0. The heat bath singles 
out a reference frame, such that for T > Galilean sym- 
metry no longer holds. In our formalism, the thermal 
contribution ut appears due to modifications of the flow 
equations for T ^ 0. We start for high k with the same 
initial values as for T = 0. As long as fc 3> ttT the flow 
equations receive only minor modifications. For k « ttT 
or smaller, however, the discreteness of the Matsubara 
sum has important effects. We plot in fig. [10] the density 
as a function of T for fixed a — 1. 

In fig. [TT]we show n{a), similar to fig. [5l but now for 
different a and T. For T — the scattering length sets 
the only scale besides n and cr, such that by dimensional 
arguments a^a — f{a'^n). Bogoliubov theory predicts 

fix) = 8T:xil + ^x^/^). (73) 




FIG. 11: (Color online) Density n for different temperatures 
and scattering length. We plot n{a) in arbitrary units, with 
A = 10'^, and for a scattering length a = 2 • 10~* (solid and 
dotted), a — 10"* (dashed and dashed-dotted). The temper- 
ature is T = (solid and dashed) and T = 1 (dotted and 
dashed-dotted). 



The first term on the r.h.s. gives the contribution of the 
ground state, while the second term is added by fluctu- 
ation effects. For small scattering length a, the ground 
state contribution dominates. We have then a a for 
n — I and a/n can be treated as a small quantity. For 
T 7^ and small a one expects a = g{T/n'^/^)an. The 
curves in fig. [TT] for T = 1 show that the density, as a 
function of cr, is below the curve obtained at T = 0. This 
is reasonable, since the statistical fluctuations now drive 
the order parameter po to zero. At very small a, the 
flow enters the symmetric phase. The density is always 
positive, but for simplicity, we show the density as a func- 
tion of (7 in figure [TT] onlv in those cases, where the flow 
remains in the phase with spontaneous U{1) symmetry 
breaking. 

For temperatures above the critical temperature, the 
order parameter po vanishes at the macroscopic scale and 
so does the condensate density ric = po = ■j^Po- The 
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FIG. 12: Order parameter po/n as a function of the dimen- 
sionless temperature T/{n^^^) for scattering length a = 10~*. 
Here, we varied T keeping a fixed. Numerically, this is equiv- 
alent to varying a with fixed T. 



density is now given by a thermal distribution of particles 
with nonzero momenta. Up to small corrections from the 
interaction aT, it is described by a free Bose gas, 



2^3/2 



(47r)3/2 



53/2(e''^), 



with the "Bose function" 
1 



9p{z) 



r(p) Jo 



(74) 



(75) 



In figure [T^ we show the dimcnsionless order param- 
eter po/n as a function of the dimcnsionless tempera- 
ture T/n^/'^. This plot shows the second order phase 
transition from the phase with spontaneous U{1) sym- 
metry breaking at small temperatures to the symmetric 
phase at higher temperatures. The critical temperature 
Tc is determined as the temperature, where the order 
parameter just vanishes - it is investigated in the next 
section. Since we find (^ — 1) <C 1, the condensate frac- 
tion ric/n = Po/n — po/{An) as a function of Tjr?!'^ 
resembles the order parameter po/n. We plot ^ as a 
function of T jr?!"^ in fig. 1131 Except for a narrow region 
around Tc, the deviations from one remain indeed small. 
Near the gradient coefficient A diverges according to 
the anomalous dimension, A ^ ^ with r\ the anomalous 
dimension. The correlation length f diverges with the 
critical exponent v ^ ^ ^ \r — Tc\^'^, such that 



A 



\T~Tr, 



-7]iy 



(76) 



Here, 77 and v are the critical exponents for the Wilson 
Fisher fixed point of the classical three-dimensional 0{2) 
model, 77 = 0.0380(4), 1/ = 0.67155(27) |i,[l3,[l3- 

In figure [TJ] we plot po jn as a function of the chemi- 
cal potential a for different temperatures and scattering 
lengths. We find, that po/n= 1 is indeed approached in 
the limit T — > 0, as required by Galilean invariance. All 
figures of this section are for u = 0. The modifications 
for u 7^ are mainly quantitative, not qualitative. 
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FIG. 13: Order parameter divided by the condensate density 
A — po/ric, as a function of the dimensionless temperature 
T /[n^^^), and for scattering length a = 10~*. Here, we varied 
T keeping a fixed. Numerically, this is equivalent to varying o 
with fixed T. The plot covers only the phase with spontaneous 
symmetry breaking. For higher temperatures, the symmetric 
phase has po = Uc = 0. The divergence of A for T ^ Tc 
refiects the anomalous dimension ry of the Wilson-Fisher fixed 
point. 
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FIG. 14: (Color online) Order parameter divided by the den- 
sity, po/n, as a function of the chemical potential. We use 
arbitrary units with A = 10'^. The curves are given for a 
scattering length a = 2 ■ 10~^ (solid and dotted), a = 10"'* 
(dashed and dashed-dotted) and temperature T — 0.1 (solid 
and dashed) and T — 1 (dotted and dashed-dotted). At zero 
temperature, Galilean invariance implies po = which we 
find within our numerical resolution. 



VII. CRITICAL TEMPERATURE 

The critical temperature Tc for the phase transition be- 
tween the superfluid phase at low temperature and the 
disordered or symmetric phase at high temperature de- 
pends on the scattering length a. By dimensional rea- 
soning, the temperature shift ATc — Tc{a) — Tc{a — 0) 
obeys ATc/Tc ~ an^^^. The proportionality coefficient 
cannot be computed in perturbation theory [isj. It de- 
pends on V and we concentrate here on u = 0. Monte- 
Carlo simulations in the high temperature limit, where 
only the lowest Matsubara frequency is included, yield 
ATc/Tc = 1.3 an^/^ Within the same setting, renor- 
malization group studies [20I [2l| yield a similar result. 
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FIG. 15: Dimensionless critical temperature Tc/{'n?^^) as i 
function of the dimensionless scattering length an^^^ (points) 
We also plot the linear fit IST^jTc = 2.\a-n}l'^ (solid line). 



for composite bosons see [22|. Near Tc^ the long wave- 
length modes with momenta ]? <C (ttT)^ dominate the 
" long distance quantities" . Then a description in terms 
of a classical three dimensional system becomes valid. 
This "dimensional reduction" is achieved by "integrating 
out" the nonzero Matsubara frequencies. However, both 
ATc/Tc and n are dominated by modes with momenta 
•p^ Ri [t^TcY such that corrections to the classical result 
may be expected. 

We have computed Tc numerically by monitoring the 
zero of po, as shown in fig. [T^ pq{T Tc) 0. Our 
result is plotted in fig. [T21 In the limit a — + we find 
for the dimensionless critical temperature Tc/(n^/^) = 
6.6248, which is in good agreement with the expected 
result for the free theory Tc/{n'^''^) = .,,1^,,,, = 6.6250. 



C(3/2)^/3 

For the shift in Tc due to the finite interaction strength 
we obtain 



AT, 



K an 



1/3 



2.1. 



(77) 



We expect that the result for k depends on the trunca- 
tion and may change somewhat if additional higher order 
couplings are included. 



VIII. SOUND VELOCITY 

The macroscopic sound velocity vs is a crucial quantity 
for the hydrodynamics of the gas or liquid. It is accessi- 
ble to experiment. As a thermodynamic observable, the 
adiabatic sound velocity is defined as 



1 dp 
Mdn 



(78) 



where M is the particle mass (in our units 1/M = 2), 
p is the pressure, n is the particle density, and s is the 
entropy per particle. It is related to the isothermal sound 



velocity vt by 



4 



1 

M 



(dp 




dp 


dT 






H 

T 


" df 







, dp 

'dT 



dT 
, dn 



where we use our units 2M 
of state" p{T, n) and 

s(r, n) 



(79) 

1. One needs the " equation 



N 



ldp_ 
ndf 



By dimensional analysis, one has 

T 



p = n^'^T{t,c), 



t = 



c = an^/^. 



(80) 



(81) 



with c) = 47rc (in Bogofiubov theory), and T{t, 0) — 
SS^l^f5/2 ^£j.gg theory), such that for small c 



C(5/2)_^5/2 



(47r)3/2 



9it)c 



(82) 



At zero temperature the second term in eq. (j79p van- 
ishes, such that Vs = Vt- For the isothermal sound ve- 
locity one has 



,dp 

dn 



_ ^dp 




( dn 
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rp da 
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\d^ 





We can now use 
dp I 

and infer 



da 



'dcrU{po) = n 



Vrp 



dlnr. 
~d^ 



(83) 



(84) 



(85) 



One may also define a microscopic sound velocity cg, 
which characterizes the propagation of (quasi-) particles. 
At zero temperature, where we can perform the analytic 
continuation to real time, we can calculate the micro- 
scopic sound velocity from the dispersion relation uj{p) 
(with p = \p\). In turn, the dispersion relation is ob- 
tained from the effective action by setting det(G~^) = 0, 
where is the full inverse propagator. We perform 
the calculation explicitly in app. [Eland find 



52 



2XpQ 



V. 



(86) 



The Bogofiubov result for the sound velocity is in our 
units 



c| — 2XpQ = Idiran. 



(87) 



In three dimensions, the decrease of S is very slow and the 
coupling V remains comparatively small even on macro- 
scopic scales, cf. fig. |31 We thus do not expect measur- 
able deviations from the Bogoliubov result for the sound 
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FIG. 16: (Color online) Dimensionless sound velocity 
Cs/{n^^^) at zero temperature, as a function of the scattering 
length an^/^. Within the plot resolution the curves obtained 
by varying a with fixed n, by varying n with fixed a, and the 
Bogoliubov result, Ca = \/T67r(an)^^^, coincide. 



velocity at T = 0. In figure (THl we plot our result over 
several orders of magnitude of the dimensionless scatter- 
ing length and, indeed, find no deviations from Bogoli- 
ubovs result. 

We finally show that for T = the macroscopic and 
microscopic sound velocities are equal, vs = vt = cs- 
For this purpose, we use 



dn I 



da da 

(88) 



From the minimum condition dpU = 0, it follows 



da 



dpdJJ_ 

diu 



a 



(89) 



Combining this with the Ward identities from app. [K[ 
namely d^U = —2VpQ and a = dpdcU — —S, valid at 
T = 0, it follows that the macroscopic sound velocity 
equals the microscopic sound velocity 



4(T = 0) 



(90) 



scattering length a and a dimensionless coupling v pa- 
rameterizing the classical dispersion relation. For T = 
and w = we find very good agreement with the Bo- 
goliubov theory. As T increases, the condensate melts 
at a critical temperature Tc, which exceeds the one for 
the free theory by ATc/Tc — 2.1an^/^. We find a sec- 
ond order phase transition in the universality class of the 
three-dimensional 0{2) model, with the associated uni- 
versal critical exponents. 

We have found an upper bound for the scattering 
length, which is of the order of the microscopic length 
scale A^^. This indicates the breakdown of the point- 
like approximation. For atom gases, A~^ typically cor- 
responds to the range of the Van der Waals interaction. 
The upper bound is a result of strong fluctuation effects 
in presence of a large microscopic interaction. Even for 
an arbitrary strong microscopic interaction, the quan- 
tum fluctuations reduce the renormalized coupling, which 
corresponds to the physical scattering length in vacuum. 
For non-relativistic atoms with pointlike interactions, a 
strong interaction regime can therefore only be realized 
at high density, n ~ A'^. Superfluid liquids, like ^He are 
good examples for such systems. First functional renor- 
malization studies for such systems have already been 
performed 23]. Our extended treatment overcomes sev- 
eral problems of this early approach and we will apply it 
to ^He in the future. Except for the vicinity of the phase 
transition, however, the detailed microphysics may play 
an important role, since all relevant length scales are of 
the order of the microphysical length scale A~^ . Other in- 
teresting extensions concern lower dimensional systems. 
We look forward to a unified functional renormalization 
description of bosons with pointlike interactions, for ar- 
bitrary dimensions and arbitrary interaction strength. 



APPENDIX A: DERIVATIVE EXPANSION AND 
WARD IDENTITIES 



IX. CONCLUSIONS 

The use of functional renormalization yields a quanti- 
tative description of the whole phase diagram for dilute 
non-relativistic bosons in three dimensions. This can de- 
scribe a gas of ultracold bosonic atoms. More gener- 
ally, our results can also be applied to quantum phase 
transitions or low temperature phase transitions, when- 
ever the most relevant excitations correspond to bosonic 
quasi-particles. For this reason we deal with a general 
dispersion relation, involving in the classical propagator 
terms linear and quadratic in the frequency. As a func- 
tion of temperature T and effective chemical potential a, 
we have computed the pressure p, the density n, conden- 
sate poj superfluid density po, and the sound velocity eg, 
in dependence on two system parameters, namely the 



We use a derivative expansion for the truncation of the 
effective average action with derivative operators up to 
four momentum dimensions 



Ffe = / {Uip,a) 



+ -Z,{p,a) {(f>*drcj)^^drr) 

+ \z2{p.a) (r(-A)<^ + (/.(- A)r) 

+ \vi{p,a) {cj>*{-dl)cj, + ^{-dlW) 
+^2 (p, o) ir (<9. A)0 - (/.(a, A)0* ) 
1 



-V^ip.a) (0*(-A2)0- 



2\ J,* 



) } (Al) 
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Here, we employ the renormalized fields 

p = = Ap = A(j)*4i 



(A2) 



and coupling functions J7, Z^, Vi. We fix the wave func- 
tion renormalization factor A such that Zi{po,cro) = 1. 
Terms of the form p(— A)p or p{—d^)p are not included 
here, since they are expected to play a sub-leading role. 
For a systematic derivative expansion they have to be 
added - the terms with up to two derivatives can be 
found in In terms of dimensions, the operator Bt 
counts as two space derivatives for the nonrelativistic 
model with 1^ = 0, while it counts as one space dimen- 
sion for the relativistic model with 3 = 0. We expand 
the fc-dependent functions U{p,a), Zi{p,a), Z2{p,a), 
Vi{p,(j),V2{p,a) and ^3(^,17) around the fc-dependcnt 
minimum po{k) of the effective potential and the k- 
independent value of the chemical potential ctq that cor- 
responds to the physical particle number density n. For 
example, with Zi — Zilpo^ao), one has 



Ziip,cr) = Zi + Z[{po,aQ){p - Pa) 
+z['^\po,cTa){a -aa). 



(A3) 



Let us concentrate on the non-relativistic model where 
S = 1, V = in the microscopic action. At zero temper- 
ature, we can perform an analytic continuation to real 
time r — it. The microscopic action ([1]) is then 



dt / (Tx 



|0*(_ja^_^_A)0 iA(0*0)}. (A4) 



In addition to space translations, rotations and time 
translations, two further symmetries constrain the possi- 
ble forms of the couplings in F. In order to derive these 
constraints, we extend eq. (|A4[) to a t-dependent source 
a{t). First, there is a semi-local U{1) symmetry of the 
form 



a ^ (J + dta. 



(A5) 



This holds since the combination {—idt — cr) acts as a 
covariant derivative. In addition, we have the invariance 
under Galilean boost transformations of the fields 

(j)*{t,x) -> = e^(«"*-'?~V*(i,^-2<7i)-(A6) 

While the invariance of the interaction term under this 
symmetry is obvious, its realization for the kinetic term is 
more involved. Performing the transformation explicitly, 
one finds 



(f>*A<j)~ 
(j)*idt(j) 



such that indeed the combination 
i9t-f A 



(A8) 



leads to this invariance. On the other hand, the validity 
of the Galilean symmetry for an effective action guaran- 
tees that only the combination (|A8[) or powers of this 
operator act on (/>. An operator of the form [idt + 7A) 
with 7^1 would break the symmetry. (Note, that Ap 
is also invariant.) 

Both the semi-local C/(l) symmetry and the Galilean 
symmetry are helpful only at zero temperature. At 
nonzero temperature, the analytic continuation to real 
time is no longer useful. An analog version of the semi- 
local C/(l) transformation for Euclidean time r would 
involve the imaginary part of the chemical potential cr, 
which has no physical meaning. The dependence of phys- 
ical quantities on ct -I- a* is not restricted. In addition, 
the Galilean symmetry is broken explicitly by the ther- 
mal heat bath. 

Combining semi-local C/(l) symmetry and Galilean 
symmetry at T = 0, we find that the derivative oper- 
ators idt, A and the chemical potential term (ct — ctq) are 
combined to powers of the operator 



i?=(-z9i-(CT-CTo)-A). 



(A9) 



In addition to powers of that operator acting on only 
spatial derivatives of terms, that are invariant under C/(l) 
transformations, like pAp, may appear. Since the sym- 
metry transformations act linearly on the fields, the full 
effective action F[(/)] is also invariant. This also holds for 
the average action Ffe[0], provided that the cutoff term 
AS'fe[(/)] is invariant. We can write the effective action as 
an expansion in the operator D 



Uoip) 



+ \z{p) ir i~idt -{a-ao)-A)^ + c.c) 
+ \v{p) (0* (-«ft - (ct - CTo) - Af4> + c.c) 



-I- . 



(AlO) 



(A7) 



Performing the Wick rotation back to Euclidean time, 
we can compare this to eq. (jAip . and find for T = the 
relations 

Zl(p,CTo) = Z2{p,(To) = Z{p), 

Fi(p, CTo) = V2(p, CTo) = Vz{p,(To) = V{p), 

z["\pa,<Jo) - 2{v{po)+PoV'{po)) , 

z("^(po,fTo) = 2y(po), (AH) 
and therefore 

a = - {z{p(^) +po^'(po)) , 
Uk = Z{pq)pq, 
P = - (2Z\po)+poZ"ipo)) . (A12) 
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We next compute the inverse propagator in a constant 
background field by expanding F^. to second order in the 
fluctuations around this background. For this purpose, 
it is convenient to decompose 



1 



(r,f) = 00 + (01 (t,:?) +i02(T,f)) . 



(A13) 



The constant condensate field (po can be chosen to be 
real without loss of generality. The fluctuating real flelds 
are the radial mode 0i and the Goldstone mode 02 , and 
P = Po + V^4'o4'i + 5^1 + 5'/'2- The truncation of the 
effective average action (|Aip reads in that basis 



+ iz2(p, a) (\/20o(-A)0i + 0i(-A)0i + 02(-A)02) 

+ V2{p,a) (iV2MdrA)(f>2 + iMdrA)c^2 - iMdrA)(f, 
+ Ivsip, a) (V20o(-A2)0i + 0i(-A2)0i + 02(-A2)<^ 



(A14) 



(2) 

The inverse propagator F^ can be inferred from an expansion in second order in 0i and 02 ■ We keep the linear order 
in (T — (Toj which will be needed for the flow equation for the density. This yields 



+ \ (Zi + Z[pQ + Z["\a - do)) (i0i9r02 - i025.0i) 

+ i (l + 2Z^po + ~ ^o)) (01 (- A)0i) 

+ ^(l + Z^'^^(d-do))(02(-A)02) 

+ \ {Vi + 2y/po + vl'\o - do)) (01 (-5^)01) 

+ {V2 + Vipo + V^"\a - do)) (i0i(a,A)02 - z02(a,A)0i) 

+ \ {y^ + 2Vipo + V^^\<j - do)) (0i(-A2)0i) 

+ ^ (^3 + Vt\<J - ^O)) (02(-A2)02) |, 



(A15) 



where we dropped the argument (po, ^o) at several places 
and used the implicit rescaling condition 2^2(po:0'o) = 1- 
In our simple truncation, we take at cr = (Tq only 

5 = Zi(/9o,Cro) + 2'{(po,0-o)po, 

V = Fi(po,ao) (A16) 

into account. We neglect the contribution of the other 
couphngs, i.e. set Z'^ = V2 ^ ^ V{ = = = 0. 
As shown above, it follows from the symmetry require- 



ments at zero temperature, that Vi 



V2 



V3 



Z'f'> = 2{V + V'po) and Z'^"' = 2V. The truncation 
V2 = Vs = therefore violates the Galilean symmetry, as 
does our choice of the cutoff term ~ Rk ■ Within our ap- 
proximation, it is consistent to set z["^ — Z^^^ = 2V at 
zero temperature. Also the deviations from this relation 
at finite temperature are neglected for simplicity in this 
paper. This yields eq. (|5D|) . 



V, 
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APPENDIX B: SYMMETRIES AND NOETHER 
CURRENTS 

In the following we discuss the role of continuous sym- 
metries of the microscopic action S[(j)]. Since all these 
symmetries are linear in the fields, the full effective ac- 
tion r[0] is also symmetric. In the case that the cut- 
off term AiS'fe[0] is chosen invariant under the symmetry 
transformation in question, this also holds for the aver- 
age action Tk[^] for finite k. From Noether's theorem it 
follows that there exists a conserved current — j) 
connected with every such symmetry. If the action is 
formulated as an integral over the imaginary time r the 
conservation equation implies for the current 



drj 



ir) 



Vj = 0. 



(Bl) 



At zero temperature, we can perform a Wick rotation to 
real time, r — > it, and eq. (|Bip takes the usual form 



(B2) 



The Noether charge C — J cPxj^*^ is conserved in time, 
i.e. j^C — 0. This holds if j falls off sufficiently fast at 
spatial infinity. At finite temperature however, the situ- 
ation is different. A simple analytic continuation to real 
time is no longer possible, since the configuration space 
is now a torus with periodicity 1/T in the r-direction. 
Instead, we can integrate eq. (jBl|) over complex time r, 
giving 

VJ^vf' dT/=/^)(0)-/^)(l/r) = 0. (B3) 

From the symmetry, it now follows that there exists a 
solenoidal vector field or three component current J = 
/r J- 

A global symmetry of an action r[0] (where F could 
be replaced by S or F^ if appropriate) can be formulated 
in its infinitesimal form as 



F[0- 



(B4) 



with e independent of x. Here s is the infinitesimal gener- 
ator of the symmetry transformation. For a local trans- 
formation, where e depends on x, e = ^{x), we can expand 



..}. 

(B5) 



The global symmetry implies that the expansion on the 
r.h.s. of eq. (jBSp starts with 5^e. Here and in the 
following it is implied that e as well as its derivatives 
are infinitesimal, i.e. we keep only terms that are lin- 
ear in e. The index /i goes over (0,1,2,3), representing 
(t, a;^, a;^, a;'^) in the real time case and (r, a;^, a;^, a;'^) for 
imaginary time. Eq. (|B5p implies for arbitrary (/'(a;) 



■5m 



= 0. 

(B6) 



Our notation is for real fields and implies a summation 
over components, if appropriate. In a complex basis one 
replaces |^es0 by ^es(j) + -^escj)* . 

Eq. (IB6|) is valid for all field configurations (f> and not 

^ = 0. In 



only for those that fulfill the field equation 
consequence, the integrand is a total derivative 



= -d^ (i^e + li^'^d^e +...). (B7) 
We can now specialize to = df^d^e = . . . = and find 

^.0(a:) = ~d,r. (B8) 

This defines the Noether current j^. For solutions of the 
field equation, = 0, the current is conserved, 

For a given x we can also specialize to 



e{x) = 0, df,e{x) ^ 0, d^d.eix) = 0, 
which leads to 



(B9) 



(BIO) 



This process can be continued, leading us to a whole 
tower of identities for the conserved current j'^. 

If the action T[(j3\ includes derivatives only up to a finite 
order n, i.e. can be written in the form 



m 



£(0, 9(/.,99(/.,...,9("V), 



(Bll) 



the expansion on the right hand side of (|B5P only con- 
tains terms up to order d'^^^e such that the tower of equa- 
tions for j'^ can be solved. Moreover, for homogeneous 



situations where 



sr[4>] 



is solved by a constant 4>, the sec- 



ond term on the r.h.s. of eq. (jB10|) vanishes since it 
includes a derivative. We have then — . 

A convenient way to find the local currents employs 
parameters e(x) that decay sufficiently fast at infinity 
such that we can partially integrate eq. (jB6|) 



This yields the local identity 



-- 0. 
(B12) 

(B13) 



An expansion of the l.h.s. in derivatives often yields sub- 
stantial information on etc. by inspection. 

Our construction yields a unique conserved local cur- 
rent 2^ for every generator of a continuous symmetry. 
We note, however, that aj^ + is also a conserved local 
current if a and are independent of x. This remark is 
important if we want to associate with the current for 
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a physical quantity. A rotation invariant setting implies 
6* = 0, but 6" and a may differ from zero. 

After these general considerations we now specialize to 
nonrelativistic real time actions of the form 



dt / d'^xCi^, {idt + A)0, {idt + A)V, • • • )■ 



(B14) 

We assume, that F invariant under the same symmetries 
as the action (|A4p . From the symmetry under time trans- 
lations 



C 



edtC 



(B15) 



we find a conserved current {jE)^ ■ Up to a possible addi- 
tive constant its i-component is the energy density, while 
the spatial components describe energy flux density. The 
multiplicative constant a gets fixed if we choose the units 
to measure energy. The choice h = \ corresponds to 
a = 1. Similarly, the invariance under spatial transla- 
tions 



(B16) 



implies a conserved current {jM)i for each spatial di- 
rection i = 1,2,3. Up to an additive constant {bM)i 
the i-component is the conserved momentum density. 
Pi — {jM)i + , while the spatial components can 

be interpreted as a momentum flux density, with the di- 
agonal components (jA/)i describing pressure. 
From the global U{1) symmetry 



C 



+ e(sc)'/' 
0* -t- e{sc)4 
C 



(f> — iecj) 



(B17) 



we can infer the conservation of the current (jc)^ as- 
sociated to the conserved particle number. In order to 
identify the total particle number with the charge of this 
current, J d^x{jc)'^ , we need to fix a possible multiplica- 
tive constant a. For this purpose, we use the Galilean 
boost invariance, described already in eq. (|A6|) . It reads 
in its infinitesimal form 

4> 4> -\- e''(sG)i4> = 4> + '2^^tdi4> ~ ie^Xi4> 

L £ + e'df,{2 5'^t£), (B18) 

and the conserved charge of (sg) is the center of mass, 
again up to an additive constant. The generator {sq) can 
be decomposed as 

{sG)^=X,{sc)-2t{sM)^. (B19) 

This implies for the current 

{jG)t = xUjcr-2t{jM)'^. (B20) 



Specializing to the t-component, identifying the momen- 
tum density pi = -I- (6M)i' and reintroducing the 
particle mass 2M = 1 we find 



M 



(B21) 



From this we can conclude that up to an additive con- 
stant (jc)° is the particle density n — {jc)° + (bc)^ . 

For the effective action (jA10[) we find for a — (Jq and 
constant (t>{x) = ^/po the current 



(jc)" = Z{po)po. 



(B22) 



Using the normalization condition Z(pq) = 1, this gives 
Uc)^ = Po- At zero temperature, this is the particle 
density and the additive constant (&c)° vanishes. At 
nonzero temperature we can compare to eq. (j72p and 
find (&c)° = riT. 

For completeness we also mention the symmetry under 
spatial rotations 



-> (l){t,R-^x) 
C{t,x) -> C{t,R-^x), 



(B23) 



with orthogonal matrix i? J- = (e"'"')^, generators {Ji)\ = 
iSijk, and Sijk the antisymmetric tensor in three dimen- 
sions. The infinitesimal transformation reads 

(/)(t, x) (j>{t, x) + ■q'erjkx^dj(l){t, x) 

C{t,x) ^ C{t,x) + r,'di{e,jkx''5]C). (B24) 

The time component of the conserved current is, 
of course, the angular momentum density. 



APPENDIX C: FLOW EQUATION FOR 
EFFECTIVE POTENTIAL 

We derive the flow equation for the effective potential 
by evaluating the flow equation for the average action 
for constant fields. Inserting a real constant field 
(f){x) ~ yfp one finds for J7 = Ffc/fi the flow at fixed p 

dtU{p,a) - rjpU' + C(p,a), 

"'0 



gi+ff2 + 2(Ui+pU/K 
h^Lol + (51 + (Fi + 2pVl)Lol){g2 + V^^l) ' 



(CI) 



Here, d is the number of spatial dimensions and we use 
the abbreviations 

gi = e -m^ + {Z2-l + 2pZ'^)p' ~{V^ + 2pVl)p^ 

+U' + 2pU", 
52 - e-m^ + {Z2-l)p'-V^v^ + U\ 
h = Z,+pZ[-{V2 + pV^)p'', 
ujn — 2T:Tn, 

Vd = (2''+V/2r(d/2))^i. (C2) 
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We dropped the arguments (p, a) at several places on the 
right hand side. Primes denote derivatives with respect 
to p. In the phase with spontaneous symmetry breaking, 
we have — dtm? = 0. 

The Matsubara sums over n can be carried out by 
virtue of the formulae 



E 



1 




tt) — \/b — dcoi\\{\l — i — it) 



Vb + dcoi\\{\j^^—TT) — Vfo — dcoi\\{\J—^—TT)j , 

(C3) 



with d = — 4ac. This brings us to 



dpp 



10 



E 



VB + D-^ - 2VB - D coth( 



/sat ' 



E 



2VbTD - ^B-D-^ ] coth( ^^^ -) 1 , 



(C4) 



where we introduced 



A = ViiVi+2pV;), 

B = +giVi+g2{Vi+2pVl), 

C = 31.92, 



D = ^Jb^ - AAC, 
E ^ gi+ g2- 



(C5) 



In our simple truncation with S = Zi + Z[ po, V = Vi, 
Z'^ = V2 = Vz = VI = = Fg' = 0, and at ct = ctq, the 
integrand in eq. (|C4[) becomes mostly independent of the 
spatial momentum. The integral can than be carried out 
and we find 



C(p,'^o)-(l 



,V2 



Vd 



d + 2' dD 



^ (^VB + D-^ - 2VB - £1^ coth( 



VB-D 



SAT 



E 



2VbTD-Vb~D-^ ] coth( ^^^ ) I, 



(C6) 



with 



A 
B 
C 
D 
E 



+ 2V{k^ + U' + pU"), 
{k^ + U' + 2pU"){k^ + U'), 



-iAC, 
U' + pU") 



(C7) 



That the momentum integral can be performed analyti- 
cally is a nice feature of the cutoff ([37|l . The limit T ^ 
is obtained by substituting the coth functions with unity. 

The flow of the effective potential contains a subtlety 
that can be seen in the limit Vi i) {i = 1, 2, 3), where 
we find 



dtU{p,a) = vpU' + 2vd 



- .92 



(2fc2 

f 91+92 .,,V9i92 
tt; — ^= coth — 
\2h^/gTg^ '2VhT 



dpp'^-^ 



The term 1 //i in the last line is not present if Vi is set to 
zero from the outset. If Zi is independent of p, this term 
is independent of p and gives only an overall shift of the 
effective potential. 



APPENDIX D: FLOW EQUATIONS FOR 
KINETIC COEFFICIENTS 



We show in this appendix our results for the flow equa- 
tion of the kinetic coefficients 5, A and V. We ne- 
glect all contributions from momentum dependent ver- 
tices. In other words, we use p-independent constants 
S — Zi + poZ[, A ^ Z2 and V — V\. In our truncation 
with Z2 = V2 = V3 = 0, and with the cutoff ([57]) . we can 
perform all momentum integrations analytically, leading 
us to 
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71 

32 Vd k^+'^X^po {k^ + fc2y) + (S'2 + 2k'^V) Xpo - 2V {S'^ + k^V + VXpo) uj,?- - 'iV^oJn^) 



5.5 = .5-(i-^)t5: 

n 



2 



d(fc4 + 2/c2Apo + + 2fc2l^ + 2FA/3o) 

A ^ d(fc4 + 2fc2Apo + (52 +2fc2y + 2FApo)w«2 + 1/2^^^4)2 



Here, d is the number of spatial dimensions, and a;„ are 
as in (|C2p . The Matsubara sums over n can be performed 
analyticaUy again by using (jCSp and derivatives thereof. 

In the hmit T 0, the Matsubara frequencies are 
continuous cj„ — + go and the sum becomes an integral 
T^Ylin ^ ^ I^oc'^^o- The expressions for ry and c^tS* 
in eq. (|Dip agree with those derived in j4j, while our 
result for dtV corrects an error in a first version of [4] (in 
equation (C.8)). 



APPENDIX E: PROPAGATOR AND 
DISPERSION RELATION 

The inverse propagator is given by the second func- 
tional derivative of the effective action 



r(2) 



s 



>i(-9) 



= G-'S{p-q), 
and we find from the truncation (jAip 

'H + 2J+{Vi + 2pVl)qi, -qoV2K 



qfi^/2K 

Here we use the abbreviations 



, H + Viq'o 



(El) 



(E2) 



J = pz;,f - pvip" + pu" 

2K = {Zi + pZ[ - 2{V2 + pV^)fy 



(E3) 



At zero temperature, we can analytically continue to real 
time qo ~f iLU, and find 



G- 



H + 2J-{Vi+ 2pV{)uj^ , ~iujV2K 



The dispersion relation is found from the on shell con- 
dition 



detG'i =0 



which yields 



+ 2HJ - 2 {H{Vi + pVl) + JVi + K)) 
+ Vi{Vi + 2pV[)u'^ = 0. 



(E5) 



(E6) 



The solutions for uj define the dispersion relation. We 
find two branches, according to 



. 2^ ^ L 

^'^^^ Vi{Vi+2pV() 



7- H{Vi + pV{) + JV1+K 



± \^{K + JVif + 2H {K{Vi + pVl) - JVpVl) 
+H\pVif\\. (E7) 



In the phase with spontaneous symmetry breaking, the 
(-f) branch of this solution is an "optical mode", while 
the (— ) branch is a sound mode. The microscopic sound 
vefocity is cs = ^l^^p. Using p = po, U' = 0, U" = X 
and Z2 = 1, we find 
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2Apn 



jzi+poz'ir IT/, 5*2 + 2XpaV ■ 



(E8) 



The "optical mode" has at vanishing spatial momentum 
the frequency 



= Q) = 



2Xpo 



jZi+p^Zj) 

Vi+2poV{ ViiVi+2poV{) 



(E9) 



.(E4) which diverges — > 00 in the limit Vi — > 0. 
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